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Abstract We investigate the early phase of the 13 February 2009 coronal mass 
ejection (CME). Observations with the twin STEREO spacecraft in quadrature 
allow us to compare for the first time in one and the same event the temporal 
evolution of coronal EUV dimmings, observed simultaneously on-disk and above- 
t he-limb. We find that these dimmings are synchronized and appear during 
the impulsive acceleration phase of the CME, with the highest EUV intensity 
drop occurring a few minutes after the maximum CME acceleration. During the 
propagation phase two confined, bipolar dimming regions, appearing near the 
footpoints of a pre-fiare sigmoid structure, show an apparent migration away 
from the site of the CME-associated fiare. Additionally, they rotate around the 
'center' of the fiare site, i.e., the configuration of the dimmings exhibits the 
same 'sheared-to-potential' evolution as the postfiare loops. We conclude that 
the motion pattern of the twin dimmings refiects not only the eruption of the fiux 
rope, but also the ensuing stretching of the overlying arcade. Finally, we find that: 
(1) the global-scale dimmings, expanding from the source region of the eruption, 
propagate with a speed similar to that of the leaving CME front; (2) the mass 
loss occurs mainly during the period of strongest CME acceleration. Two hours 
after the eruption Hinode/EIS observations show no substantial plasma outfiow, 
originating from the 'open' field twin dimming regions. 
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1. Introduction 



Coronal dimmings are regions in the solar corona that undergo an abrupt inten- 
sity drop roughly co-temporal with the launch of a coronal mass ejection (CME). 
They are observable in soft X-rays (SXR; [Sterling and Hudson , 1997) and most 
conspicuously in the extreme UV (EUV; e.g.. 



al. 



2000) both on the solar disk and above the limb. 



Zarro et al. , 1999 Thompson et| 



^ Space Research Institute, Austrian Academy of Sciences, 
SchmiedlstraBe 6, A-8042 Graz, Austria 
^ Kanzelhohe Observatory-IGAM, Institute of Physics, 
University of Graz, Universitatsplatz 5, A-8010 Graz, Austria 
emaih asv@igam.uni-graz.at| 



SOLA: mikl_v2_julll.tex; 4 October 2011; 0:26; p. 1 



Global-scale dimmings originate from the source region of the eruption and 



may expand across almost the entire solar disk (e.^., Zhukov and Auchere, 2004 



Mandrini et al. 2007). They follow the nearly spherically propagating fronts 



of global coronal waves (Moreton waves and EIT waves, Moreton and Ramsey 
1960[ [Moses et al\ |1997| [Thompson et al\ |1998[ ). On the other hand, small- 



scale twin dimmings, which typically exhibit a more dramatic intensity decrease, 
occur in confined areas near the ends of a pre-flare, S-shaped sigmoid structure 
(Sterling and Hudson, 1997). 



The generally accepted physical interpretation of coronal dimmings is that 
they reflect the strong density depletion of the inner corona due to the expansion 
or 'opening' of magnetic field lines in the wake of a CME. Two major lines of 
evidence support this view: (1) Measurements of Doppler shifts in emission lines 
formed in the lower transition region and corona revealed material outflows co- 
spatial wit h small-scale dimming regions, with outflow velocities of seve ral tens 
of kms-i ( [Harra and Sterling! |2QqT| jHarra et "a/] [2007|| Jin et a/.[|2QQ9l ). Inside 
global-scale dimmings mass outflows with velocities of ~ lOOkms"^ were found 
(Harra and Sterling, 2001]). (2) Based on the idea that due to mass conservation 
the mass accumulated in the leading edge of the CME has to be balanced by 
the density depletion in the dimming region, coronal dimmings have been used 



to estimate CME masses (e.^.. 


Harrison and Lyons 2000 Harrison et al. , 2003 


jAschwanden et al. 2009a|b 


Jin et al. 20091). In some cases, the consistency with 



CME masses based on white-light measurements was good (e.^., Harrison etal\ 



2003 Aschwanden et al., 2009a), in other cases, the CME mass based on EUV 



observations was 30-70% smaller than the white-light estimate {e.g., 



and Lyons 2000 Jin et al. , 2009 ) . Jin et al. ( 2009 ) attributed this discrepancy 



Harrison 



partly to transition region outflows, which refill the dimming regions. Never- 
theless, these studies strengthen the link between coronal dimmings and CMEs. 
Therefore, the investigation of the dimming process is considered as a powerful 
diagnostic of the early phase of a CME, especially with regard to the physics of 
CME onsets ([Harrison et al. , 2003). This phase cannot be examined using white- 
light observations, since the occulting disk of the coronagraph obscures both the 
Sun and the lower layers of the corona, where the CME is usually launched and 



accelerated (Temmer et al, 2010 2008). 



Previous studies on coronal EUV dimmings addressed their temporal and 
spatial relationship to the associated CME, the amount of the EUV intensity 



drop as well as the timescale and duration of the dimmings (e.^.. Sterling and 
Hudson! [T997| [Aschwanden et al.\ [19991 [Mcintosh et al.\ |2QQ7| jBewsher et al. 



2008). Mcintosh et al. (2007) examined both the temporal evolution of the size 
of particular dimming areas and the apparent motion of these regions. Due to 
observational limitations, however, in all of these studies a particular dimming 
event was observed either on the solar disk or above the limb, depending on the 
particular line of sight from the observing instrument to the source region of the 
eruption. Therefore, it was impossible to examine in one and the same event the 
temporal evolution of coronal dimmings from different vantage points. 

With the launch of the twin STEREO spacecraft {Solar- Terrestrial Relations 
Observatory^ Kaiser et al. , 2008 ) these observational limitations are temporarily 



over, and more importantly, the Sun was cooperative enough to produce a 
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CME/flare event when STEREO Ahead (STA) and STEREO Behind (STB) 
were in quadrature on 13 February 2009. The resulting unprecedented data 
set ahows us to investigate the coronal EUV dimming process in this event 
simultaneously both on the solar disk and above the limb. Due to the line- 
of-sight integration of the optically thin emission, the combination of on-disk 
and above- the- limb observations is particularly instructive. We compare both 
dimming observations and relate them to the kinematics of the associated CME. 
Finally, we combine EUV observations obtained from both the STEREO and 
the Hinode spacecraft (Kosugi et al. 2007 ) to examine the relationship between 



dimming locations and mass flows. 

In Section [2] we describe the observations as well as the data sets used and 
the data reduction procedures applied. The analysis and results are presented in 
Sections [3] and m In Section [5] we discuss the results. 



2. Observations and Data Reduction 

The 13 February 2009 CME/fiare event occurred in NO A A active region AR 
11012. At this time, the separation angle of the two STEREO spacecraft was 
91°. STA observed the event above the solar limb, while STB observed it near 
the center of the solar disk. The CME had a typical three-part structure, i.e., a 
bright front, a dark cavity, and a bright core. The flare (GOES class B2.4) started 
around 05:30 UT (GOES maximum at 05:47 UT). The event was associated with 
a global coronal wave observed in EUV ( Cohen et a/.[|2009| Kienreich et a/. [[2009 



Patsourakos and Vourlidas 2009 ) . 



Although in this paper we focus on the investigation of the early phase of the 



CME, it is interesting to note that the Heliospheric Imagers (HI, Eyles et al. 



2009 ) onboard STA tracked the CME's leading edge up to the distance of about 



1 AU. In add i tion, the STB/IMPACT ( [Luhmann eTal\ [2QQ8[ ) and PLASTIC 



(Galvin et al. 2008) instruments registered a weak shock on 18 February 2009, 



10:00 UT, followed by the complex magnetic field and bulk plasma signatures 
typically associated with an interplanetary coronal mass ejection (ICME). A 
possible encounter of the ICME with Venus was observed with the Venus Express 
magnetometer (Zhang et al. , 2006[ ). A detailed study on the ICME using the 
STA/HI observations along with the in situ measurements provided by STB and 



the Venus Express is underway and will be presented elsewhere (Mostl et al. 
[2011] ). 



2.1. STEREO Data 



From the STA and STB spacecraft we use EUV 171 and 195 A filtergrams. 



\et 



provided by the Extreme Ultraviolet Imager (EUVI; Wuelser et al. 2004 Howard 
2008[ ), to investigate the dimming process both on the solar disk and above 



the limb. In addition, we use the STA/EUVI 171 A filter in combination with 
STA/CORl white-light images to track the CME front, and thus, determine its 
velocity and acceleration profile. The observing cadence in the EUVI 171 A filter 
is 2.5 min for STA and 5 min for STB, while in the 195 A filter it is 10 min 
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for both STEREO spacecraft. The white-hght STA image cadence is 10 min for 
CORl and 15 min for C0R2. 



The SECCHLPREP routine in the SolarSoft software package (SSW; Freeland| 



and Handy, 1998) was used to reduce the images. Ah images were rescaled 



to Earth distance to compare observations from STA and STB. STB/EUVI 
filtergrams were differentiahy rotated to the same reference time. In addition, 
the gamma-log function, a compound of the logarithmic and gamma functions, 
was applied to the STB/EUVI filtergrams to enhance faint structures in the 
images while retaining detail in the bright features. 

2.2. Hinode Data 

To derive plasma flow velocities in the active region we use the earliest available 
EUV raster scan for this event, provided by the EUV Imaging Spectrometer 
(EIS; [Culhane et a/., 20071) onboard the Hinode spacecraft. The scan covers the 



period immediately after the decay phase of the CME- associated flare (07:25 - 
08:12 UT). The EIS instrument provides high-resolution EUV spectra within two 
wavebands: 166 - 211 A and 246 - 292 A, covering emission lines in the transition 
region and corona. Four slit/slot positions (1", 2", 40", and 266") allow a wide 
choice of operation modes. EIS raster scans are carried out in the following way. 
In each emission line the intensity is measured by sampling the line across the 
1" or 2" width of the slit. This sampling is carried out at up to 512 positions 
across the length of the slit. After this, the slit is moved to the next position 
of the region of interest until the entire region is scanned. Observations can be 
carried out in up to 25 lines simultaneously. Images in a particular emission 
line are constructed as follows. To obtain a pixel column of the image, at each 
position across the length of the slit, the intensity is integrated over the spectral 
line width. Repeating this procedure for all slit positions within the region of 
interest yields the entire image. 

After the decay phase of the B2.4 flare, the EIS instrument scanned the active 
region from west to east in 90 steps, each step 2" wide, with a slit length of 
320", yielding a field of view (FOV) of 180" x 320". The observations comprise 
eight spectral windows, each 32 pixels wide (i.e., the intensity in each line was 
sampled at 32 positions across the line). We use the Fe xil 195.12 A window. We 
calibrated the EIS data using the IDL routine EIS_PREP, available in the EIS 
part of the SolarSoft tree. This routine corrects for fiat field, dark current, and 
cosmic rays, and flags saturated, hot, and warm pixels. After this, we derived 
a Dopplergram from the Fe xil 195.12 A intensity image as follows. Using the 
routine EIS_AUTO_FIT, we fitted a single Gaussian to the calibrated EIS data 
to compute a velocity array from the measured line centroid positions. We then 
corrected the obtained velocity and centroid vectors for the tilt of the EIS 
slit (EIS_TILT .CORRECTION) and the orbital variation of the line centroid 
(EIS_ORBIT_SPLINE). 

2.3. Co-alignment of STEREO/EUVI and Hinode/ElS Images 

As aforementioned, EIS images, derived from raster scans, are not snapshots 
taken at a particular moment in time. When we look at such an image we look at 
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stringed-together pixel columns taken at different instants in time. Therefore, if 
the scanned region is a flaring region, the features in the FOV may have changed 
dramatically during the scan. Accordingly, the co-alignment of EIS images and 
EUVI filtergrams is a tricky issue. Even though there will be similar features 
in EIS and EUVI images taken in the same passband, the images will not look 
exactly the same. To solve the problem of co- alignment, we used a method similar 
to that described by Jin et al. (2009). We compared EIS and EUVI images at 
195 A. We set the mid-time of the raster scan as the time of the EIS image, chose 
an EUVI image observed closest to this time, visually compared it with the EIS 
image by identifying similar features, and then shifted the values of XCEN and 
YCEN in the EUVI image to roughly co-align it with the EIS raster image. 
After this, we cross-correlated both images to refine the visual co-alignment. We 
estimate the uncertainty of the co-alignment achieved with this method to ±3". 



3. Analysis 

STEREO/EUVI 171 and 195 A filtergrams reflect the emitting state of the 



coronal plasma at approx. 1 and 1.5 MK, respectively (Wuelser et al. 2004 
[Howard et a/. [[2008 ). To investigate the coronal dimming process, we calculate 



intensity light curves from these filtergrams by (1) summing up at each time all 
pixel intensities within a region of interest (ROI) and (2) dividing this sum by 
the number of pixels in this region. For a global light curve the ROI is a large 
rectangular box surrounding the dimming regions, for a local light curve the 
ROI is determined as follows. To enhance the contrast of the dimming regions 
we use base difference images, i.e., we subtract a pre-event image taken around 
13-Feb-2009 04:45 UT from all successive frames. We then count among the 
dim pixels all pixels with intensities in the range [/min, / • ^min], where /min is 
always negative, and is the intensity of the darkest pixel in a difference image 
and /e[0, 1]. For instance, for /min = —500 (in units of normalized detected 
photons) and a factor of / = 0.8, pixels with intensities in the range [—500, 
—400] are counted among the dim pixels (z.e., an 80%-threshold is used). Hence, 
the number of pixels forming the ROI of a local light curve varies in time. This 
offers the opportunity to track the size evolution of local dimmings as well. We 
normalize all light curves to the pre-event intensity level. 

The drop in the EUV intensity / during coronal dimmings is usually inter- 
preted as mass loss in the corona in the wake of a CME. For instance, jJin et 



al. (2009) showed that the EUV intensity variation in coronal dimming regions 
reflects the variation of the mass density. Therefore, we consider the intensity 
change rate dl/dt as a proxy for the mass loss rate, i.e., the time evolution of 
the mass loss. We calculate dl/dt as the time derivative of a spline fit to the 
normalized dimming light curve. 

We derive the temporal evolution of the CME velocity and acceleration (i.e., 
the profiles vcME{t) and acME(^), respectively) as follows. We (1) construct 
running difference images from STA/EUVI 171 A and CORl white light frames 
to enhance the contrast of the CME front, (2) track the fastest-moving part of 
the front, (3) fit the obtained height-time measurements by a spline fit, and (4) 
calculate the first and second time derivative of the fit curve. 
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4. Results 



4.1. On-Disk EUV Dimmings and Associated CME/Flare Evolution 
4-1-1- Flare Evolution and Duration of the Dimmings 

Figure [l] shows the source region of the CME observed with STB/EUVI 171 A. 
The CME-associated flare occurs at the site of an S-shaped sigmoid structure 
(Figure [TJ a)). During the impulsive phase of the flare, two elongated, roughly 
north-south aligned flare ribbons separate (Figurejljb)). The first postflare loops, 
resulting from the reconnection of magnetic field lines, appear at the beginning 
of the decay phase (Figure [TJc)). Furthermore, two small-scale dimming regions, 
designated as DRl and DR2, appear near the ends of the pre-flare sigmoid. 
These regions grow and exhibit an apparent motion during the decay phase 
(Figure [TJd)). The postflare loops in between the dimmings form an east- west 
aligned arcade. 

In the top row of Figure [2] we plot the STB/EUVI 195 A intensity light curve 
calculated from the FOV of the images presented in the middle row (global light 
curve). At the beginning of the presented time interval we use the full image 
cadence available (10 min) to capture the abrupt intensity drop. After 13-Feb- 
2009 08:25 UT the cadence is reduced to 40 min, as a lower cadence is sufficient 
to capture the subsequent gradual changes in the EUV intensity. The intensity 
drops down to roughly 83% of the pre-flare EUV level in approx. 30 min. The 
recovery of the EUV intensity back to the pre-event level in the event under 
study takes approx. 16 h. We note that in extreme cases, recovery times of 



coronal dimmings up to 48 h have been observed (Attrill et al. 2008). The 
middle and bottom rows show four representative EUVI 195 A direct and base 
difference images, respectively, illustrating the onset of the dimmings and the 
ensuing gradual intensity increase, as the rarefied regions are slowly refilled. 



4.1.2. Apparent Motion and Size of the Small-Scale Dimming Regions 

We use the STB/EUVI 195 A filter grams to track the apparent motion of the 
localized dimming regions DRl and DR2 {of Figure [l] panels (c) and (d)) and 
to investigate the evolution of their size. We determine the pixel size of DRl and 
DR2 in each base difference image, using a factor of / = 0.8. This 80%-threshold 
proved to be most suitable for tracing the dim areas in an image. After this, we 
calculate the 'center of gravity' (CoG) of DRl and DR2 at each time (i.e., for 
each frame we weight the location of dim pixels by their intensity, and then we 
divide the sum of these products by the number of dim pixels). Figure |3] shows 
the locations of the CoGs for both dimming regions at each time superimposed 
on a base difference image, taken during the decay phase of the flare. The CoG 
of DRl moves to the north-east, while the CoG of DR2 is directed to the south- 
west. The straight lines connecting the CoGs of both dimming regions at the 
beginning and end of the analyzed time interval demonstrate that the CoGs 
rotate clockwise around the 'center' of the flare site. Also, taking the east-west 
alignment of the postflare loops around 07:45 UT as a reference, the lines indicate 
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Figure 1. Evolution of the CME-associated flare and the small-scale coronal dimmings in 
STB/EUVI 171 A. (a) Pre-flare S-shaped sigmoid structure, (b) Bright, separating flare rib- 
bons shortly after flare onset, (c) The flrst postflare loops, resulting from the reconnection of 
magnetic fleld lines, appear at the beginning of the decay phase of the flare. The white arrows 
point to two small-scale dimming regions, designated as DRl and DR2, which appear near 
the ends of the pre-flare sigmoid structure, (d) During the decay phase, the dimming regions 
exhibit an apparent motion. The postflare loops in between the dimmings form an east-west 
aligned arcade. The dashed box shows the fleld of view (FOV) of Hinode /KIS. - North is up, 
west is to the right. 

that the observed rotation reflects the evolution of the dimmings from a 'sheared' 
to a more 'potential' configuration, as the event progresses. The successively 
appearing postflare loops exhibit the same sheared-to-potential evolution (c/. 
Figure [l] panels (c) and (d)). 

In Figure [4] we consider the apparent motion of the CoGs of the two dimming 
regions separately for the east-west (x) and north-south (y) direction. In addi- 
tion, we compare the motion of the CoGs with both the kinematics of the CME 
and the temporal evolution of the soft X-ray emission of the associated flare. The 
CoGs of the dimming regions separate in the x-direction (Figure [4]^ a)), while 
they approach each other in the ^/-direction (Figure |4jb)). Given the east- west 
alignment of the latest appearing postflare loops, the approaching in y reflects 
the aforementioned evolution of the dimming regions from a 'sheared' to a more 
'potential' conflguration. The separation in x indicates that the extent of the 
region, affected by the expelled CME, grows, as the event progresses. The most 
conspicuous movement of the CoGs is detected between 05:50 - 06:30 UT {i.e., 
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Figure 2. Evolution of the on-disk EUV dimmings. Top: STB/EUVI 195 A light curve. 
The intensity values are normalized to the pre-flare EUV flux. The shaded area highlights the 
period of the dimmings and the EUV flux recovery back to the pre-event level. The times of the 
images are marked with triangle symbols. The times of the four selected images shown below 
are marked with dashed hues. Middle: STB/EUVI 195 A images covering a 600"x500" FOV 
centered around the dimmings. Bottom: Corresponding base difference images. At 06:05:49 UT 
two distinct small-scale dimming regions are visible. They slowly disappear during the next 
16 h. 



after the impulsive acceleration phase of the CME). In this period, the CME is 
in its propagation phase, moving with a constant speed of around 350kms~^ 
(Figure [4]^ c)). After approx. 07:00 UT, the CoGs of both dimming regions seem 
to stall in place. 

In Figure |4]^d) we present the CME acceleration profile acME- The impulsive 
acceleration phase of the CME lasts for approx. 40 min and reaches its maximum 
of 250 ms~^ at 05:36 UT. Around this time, the dimming regions appear and 
the CoGs start to move. In Figure ^d) we also investigate the relationship 
between the acceleration phase of the CME and the energy release phase of the 
flare. The maximum CME acceleration occurs approx. 11 min before the flare- 
related maximum of the soft X-ray (SXR) emission. In other words, the highest 
CME acceleration occurs in the impulsive phase of the flare, during which the 
flare energy is released and particles are accelerated to nonthermal energies. 
Since observations of nonthermal emission (hard X-rays, microwaves produced 
by these particles) are not available for the event under study, we assume that 
the Neupert effect is valid in this event, i.e., we assume that peaks in the time 
derivative of the SXR flux coincide with peaks in the nonthermal hard X-ray 
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Figure 3. STB EUVI 195 A base difference image of the flaring and localized dimming regions. 
Bright east-west aligned postflare loops are visible in the center of the image. Gray contours 
trace the total dimming area, i.e., the area that exhibited a drop in intensity within the time 
interval indicated by the color bar. The contour level represents the 80%-threshold used to 
identify the dimming regions. The color-coded diamonds mark the positions of the 'center 
of gravity' (CoG) of dimming regions DRl and DR2 at each time. The dotted/dashed line 
connects the GoGs of both dimming regions at the beginning/end of the time interval. 



(HXR) emission {e.g., Neupert 1968 Dennis and Zarro 1993 Dennis et al. 



2QQ3[ [Veronig et al. , 2005). Consequently, we fit a spline curve to the GOES 



SXR light curve and take the time derivative of the fit to obtain a proxy for 
the HXR emission, and thus, a proxy for the energy release rate in the flare 
{of. dashed profile in Figure |4]^d)). The maxima in the SXR time derivative and 
CME acceleration profile are roughly co-temporal. The time difference between 
the peaks is less than 4 min. This indicates that the CME acceleration phase is 
closely synchronized with the energy release phase of the associated flare {e.g., 



Zhang et al. 2001 Maricic et al. 2007 Vrsnak et al. 2007 Temmer et al. 2010). 



In the upper panel of Figure 5] we present the temporal evolution of the size 
of DRl and DR2. Both dimming regions grow steadily until approx. 07:20 UT, 
but DR2 grows faster than DRl. After this, the size of DRl remains almost 
constant within the time interval under study, while DR2 starts to contract as 
it is slowly refilled. The bottom panel of Figure |5] shows the time history of 
the intensity change of DRl and DR2, averaged over all pixels in the respective 
dimming region (local light curves). Both regions undergo an abrupt intensity 
drop on a timescale of approx. 25-30 min. The intensity drops down to roughly 
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Figure 4. Apparent motion of dimming regions DRl and DR2 in east-west (x) and north- 
south (y) direction related to the CME kinematics and the soft X-ray (SXR) emission of the 
flare, (a) and (b): The x- and y-distance of the 'centers of gravity' (CoG) of DRl and DR2 
observed in STB/EUVI 195 A. (c): Height-time measurements of the CME front observed in 
STA/EUVI 171 A and CORl white hght plus sphne flt (rcME) and derived velocity f;cME- (d): 
Derived CME acceleration acME (solid thick), GOES12 1 - 8 A SXR flux (dotted) plus spline 
fit (-FsxR, solid thin), and time derivative of the flt (dashed). The time derivative is scaled in 
such a way that its amplitude is equal to the peak of acME- The vertical bar marks the time 
difference of 3.5 min between the maxima of acME and GOES derivative. The horizontal bar 
shows the period of the EIS raster scan. 
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Figure 5. Temporal evolution of the dimming areas. Top: Size of dimming areas DRl and 
DR2 observed in STB/EUVI 195 A. Bottom: Normalized intensity, averaged over all pixels in 
the dimming region. The horizontal bar shows the period of the EIS raster scan. 



20% of the pre-flare EUV flux and remains virtually unchanged during the next 
2.5 h. 

4- 1-3. Dimming Process and Mass Loss Rate Compared to the CME/Flare 
Progression 

In the top panel of Figure |6] we plot: (1) the SXR emission profile of the flare 
along with its time derivative; (2) the normalized STB/EUVI 171 and 195 A 
intensity profiles, calculated over a FOV of 640" x 596" surrounding both the flare 



SOLA: mikl_v2_julll.tex; 4 October 2011; 0:26; p. 11 



site and the small-scale twin dimmings (global light curve); and (3) the spline 
fits / to the data points. During the impulsive phase of the flare, the intensity 
drops down to a level of approx. 90% of the pre-flare EUV flux in 30 min. In 
the bottom panel of Figure [6] we investigate the temporal relation between the 
acceleration phase of the CME (represented by the acceleration profile acME) 
and the EUV intensity change rate dl/dt^ which is considered as a proxy for 
the evolution of the mass loss rate {cf. Section [3]). In both wavelengths, the 
intensity change is most conspicuous during the acceleration phase of the CME. 
The peak of acME appears at 05:36 UT. The peaks in the mass loss rate occur a 
few minutes later, namely, at 05:41 UT in the 195 A filter and 05:43 UT in the 
171 A filter, respectively. This indicates that large amounts of mass are removed 
from the analyzed region while the CME is expelled from the Sun. Moreover, 
the timing of the highest mass loss coincides roughly with the maximum of the 



SXR time derivative (05:40 UT, cf. Section 4.1.2), suggesting that the energy 



release in the flare is synchronized with the mass loss due to the CME. 
4.2. Plasma Flows in the CME Source Region 

Large amounts of mass expelled from the Sun in the course of the CME should be 
reflected in upflowing plasma at the source region of the eruption, in particular, 
at the locations of the on-disk dimming regions DRl and DR2, observed with 
STB/EUVI 195 A {cf. Figure [3|. To measure the plasma flow velocities in the 
dimming regions we use the first available Hinode /EIS raster scan and derive a 
Dopplergram from observations in the Fe xil 195.12 A emission line. 

Figure[7|a) shows an STB/EUVI 195 A intensity image taken around the mid- 
time of the EIS raster scan. The FOV is the same as in the EIS instrument {cf. 
dashed box in Figure [TJd)). The contours mark the total area of DRl and DR2 
(white) and the EIS flow velocities between —25 and +15kms~^ (blue/red). 
The same contours are superimposed on the Dopplergram in Figure ^h) . In 
DRl we detect only downflows, and these downflows are relatively weak (around 
5kms~^). The strongest downflows with velocities around 15 kms~^ are ob- 
served in the arcade of postflare loops, predominantly at the loop footpoints 
{cf. red contours in Figure [7|^ a)). In DR2 we find both downflows and upflows, 
but the upflow velocities are low (less than approx. 5kms~^). The strongest 
upflows with velocities around 20-25 kms~^ are co-spatial with a large-scale, 
coronal loop structure south-east of DR2 {cf. dark-blue contours in Figure [t]^ a)). 
This implies that during the EIS raster scan {cf. Figure ^d)) the strongest 
upflows of the 1.5 MK plasma seem to be rather flows along large, closed coronal 
loops than outflows originating from 'open' field regions. This is consistent with 
the finding presented in Section |4.1.3 The main mass loss, inferred from the 



STEREO/EUVI intensity drop occurs between roughly 05:25 and 05:55 UT, 
i.e., approx. 90 min before the beginning of the EIS raster scan {cf. Figure [6] 
and [5|. During the time interval of the EIS observations, the mass loss rate 
approaches zero, indicating that the coronal mass loss in the wake of the CME 
has stopped. 
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05:00 05:20 05:40 06:00 06:20 06:40 
Stort Time (l3-Feb-09 04:45:49) 

Figure 6. Top: GOES12 1 - 8 A soft X-ray flux (SXR, solid thick), GOES time derivative 
scaled to the SXR peak flux (solid thin), normahzed STB/EUVI 171 and 195 A hght curves 
integrated over the FOV of Figure [oj), and spline flts /in and /195 to the data points. The 
thick X -symbol around 05:32 UT highlights the temporary increase in the 171 A intensity, 
caused by the onset of the flare. The vertical line at 05:39:30 UT indicates the maximum of 
the GOES derivative. Bottom: CME acceleration proflle acME (solid) and intensity change 
rates (time derivatives of the spline curves /171 and /195). The vertical lines at 05:40:49 and 
05:43:19 UT indicate the times of the highest change in intensity for both fllters. The vertical 
bar marks the time difference of 1.3 min between the maximum of the GOES derivative and 
the minimum of dliQ^/dt. 

4.3. EUV Dimmings Observed Above the Solar Limb 

We use the STA/EUVI 195 A filtergrams to investigate the dimming process 
above the solar limb. Figure [sj^a) shows the extent of the dimming region, 
designated as DR-L, at 06:05 UT (i.e., in the early decay phase of the flare). 
A 70%-threshold was used to determine the pixel size of DR-L at each time 
(/ = 0.7). In Figure (sjb) we present both a local and a global normalized 
intensity light curve. The global profile is calculated from the entire FOV of 
Figure [s]^ a), while the local one is derived from all pixels forming DR-L at each 
time. The spline fits to the profiles show that the EUV intensity decrease is more 
distinct inside DR-L than in the entire FOV (intensity drop to 32% vs. 86% of 
the pre-CME flux in 25 vs. 40 min). After 06:00 UT, a slight gradual intensity 
increase is detected in both profiles. Figure (sjb) also shows the size evolution 
of DR-L. The dim area grows steadily during the period of intensity drop and 
reaches its maximum size around 06:05 UT. After this, the area contracts slowly 
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Figure 7. (a): STB/EUVI 195 A intensity image taken at 07:45:49 UT (i.e., around the 
mid-time of the Hinode /FjIS raster scan), showing the region of the EIS FOV (cf. dashed 
box in Figure [ijd)). The dashed hue marks the position of the EIS sht at 07:49:11 UT. (b) 
Dopplergram derived from observations in the Hinode /FjIS Fe xil 195.12 A emission hue. 
Positive and negative velocities correspond to red and blueshifts, respectively. The time axis 
at the top gives the beginning, mid-time, and end of the EIS raster scan. - Both panels: 
Black/white contours mark the total area of the STB/EUVI 195 A dimming regions DRl 
and DR2 {cf. Figure [sj. The contour level represents the 80%-threshold used to identify the 
dimming regions. The colored contours mark Hinode /EIS 195.12 A plasma flow velocities with 
contour levels of [-25, -20, -15, -10 ] kms'^ (blue) and [+10, +15] kms"! (red). 



for about an hour, and then its size remains virtuahy unchanged within the time 
interval under study (i.e., until 08:25 UT). 

Figure [sjc) shows the time derivative of the local intensity profile /dr-l- The 
curve progression indicates that the intensity drop, and therefore mass loss, is 
most distinct during the period highlighted by the shaded box in Figure [sjc). 
The width of this box is the FWHM of the CME acceleration profile shown 
in Figure ^d) and in the bottom panel of Figure [gJ We use the FWHM to 
represent the period of strongest CME acceleration. The end of this period 
roughly coincides with the end of the impulsive phase of the flare {cf. vertical 
solid line at around 05:48 UT in Figure [sj^c), marking the GOES maximum). 
The peak time of the intensity change rate at 05:37 UT is close to the middle of 
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Figure 8. EUV dimmings above the limb, (a): STA/EUVI 195 A difference image with a 
pre-CME frame around 04:45 UT subtracted from a post-eruption frame at 06:05:30 UT. The 
sohd contour enclosing DR-L represents the 70%-threshold used to identify the dimming region. 
The dashed arc marks the solar limb. The solar disk has been artificially occulted, (b) Solid: 
temporal evolution of the size of DR-L observed in STA/EUVI 195 A. Dotted: spline fit /fov 
to the A-symbols, which represent the normalized intensity calculated from the FOV of panel 
(a) at each time. Dashed: spline fit /dr-l to the 0-symbols, which represent the normalized 
intensity calculated from all pixels forming DR-L at each time, (c): Time derivative of the 
spline fit /dr-l- The gray vertical bar highlights the period of strongest CME acceleration, 
represented by the FWHM of the CME acceleration profile shown in Figure ^d). The vertical 
lines mark the time of the minimum of the intensity change rate at 05:37:30 UT (dashed) and 
the time of the GOES maximum of the fiare at 05:47:39 UT (solid). The horizontal black bar 
highlights the period of the EIS raster scan. 



the FWHM box (05:36 UT). After the period of strongest CME acceleration the 
intensity change rate approaches zero, in agreement with the intensity change 
rate derived from the on-disk observations {cf. Section 4.1.3). This indicates that 
the mass loss is synchronized with the CME acceleration phase and substantial 
only during the period of strongest CME acceleration. 



4.4. CME Velocity in the Initial Stage of the Eruption 

In Figure Qc) we presented the evolution of the CME velocity, derived from 
combined STA/EUVI and CORl above- the- limb observations, extending over 
the period between roughly 05:20 and 06:45 UT. In the following, we exploit the 
availability of both above- the- limb and on-disk EUVI observations in one and 
the same event to compare two methods of determining the velocity of the CME 
in the initial stage of the eruption (i.e., until approx. 05:56 UT). 

In the first method we use above-the-limb observations. We apply the pro- 
cedure described in Section |3j i.e., we track the fastest part of the CME front 
m STA/EUVI 171 A running difference images and calculate the velocity from 
a linear fit to the height-time measurements. Figure [9ja) shows an example of 
such a measurement. To estimate the error we repeat the measuring procedure 
several times and determine the mean and standard deviation of the calculated 
velocities. Thus, we obtain a linear CME velocity of 212 ± lOkms"^ between 
05:26 and 05:46 UT. 
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Figure 9. (a): STA/EUVI 171 A running difference image. Dashed line: direction, in which 
the CME front was tracked; x-symboh one of the height-time measurements; sohd arc: solar 
limb, (b): STB/EUVI 171 A filtergram taken at the same time as (a). The dashed line from 
the flare site to the most remote vertex of the FOV marks the distance used to estimate the 
average CME velocity. Arrows point to dimming regions DRl and DR2. FOV: 640" x 596". 
(c): Solid/dotted: STB/EUVI 171/195 A flux integrated over the FOV of panel (b) and 
normalized to the pre-flare intensity. The A-symbols represent the averaged STA height-time 
measurements of the CME front. The dashed line is the linear flt to the data points inside the 
gray vertical bar that highlights the period At of the EUV intensity drop. 



In the second method we use STB/EUVI 171 and 195 A on- disk observations 
of the region surrounding both the flare site and the smah-scale twin dimmings. 
To select an appropriate FOV, i.e., a FOV covering the range over which the 
outermost border of the expanding global-scale dimmings is discernible, we vi- 
sually track these dimmings in difference images. Figure [9]^b) shows the FOV 
selected according to this criterion. We associate the expanding global-scale 
dimmings with the 1 MK and 1.5 MK plasma, rarefied or ejected in the course 
of the CME (c/., Aschwanden et at. , 1999). We calculate the intensity profiles for 
both wavelengths, integrating over the entire FOV of Figure [gfb) (global light 
curve) . We then normalize both profiles to the pre-flare intensity. Figure |9jc) 
shows that the intensity drops down to a level of 87-90% of the pre-flare EUV 
flux on a timescale of At ^ 30 min (05:26 - 05:56 UT). This is the period 
over which the CME clears the FOV, because afterwards the intensity level in 
both wavelengths stabilizes. The average velocity of the CME in this period is 
calculated as vqme = dj ISi^ where d (~ 345,000 km) is the distance shown in 
Figure [9]^b) . To estimate At we use both the EUVI 171 and 195 A light curves 
and determine the times when the intensity drop in both profiles starts and ends 
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{cf. shaded box in Figure [9]^c) , highhghting the time interval At). This yields 
an average speed of ^cme ^ 192kms~^. Due to the image cadence of 5 min for 
STB 171 A and 10 min for STB 195 A, however, At is likely to be overestimated, 
since the actual intensity drop may start later, namely, between the time of the 
last pre-event intensity image and the time of the first 'intensity-drop' image. 
Similarly, the actual intensity drop may end earlier, namely, between the time 
of the last 'intensity-drop' image and the time of the first image with stabilized 
intensity. Assuming that the overestimation of At is up to 10 min, we obtain an 
average speed of ^cme ^ 265kms~^. Taking the mean of both velocity values 
we obtain an average speed of '^cme ^ 228 ± 36kms~^. 

Both the 'above- the- limb' method (CME tracking) and the 'on-disk' method 
(expanding global-scale dimming) yield comparable results. 



5. Discussion 



Taking advantage of the quadrature of the STEREO Ahead and Behind space- 
craft on 13 February 2009, we presented the first direct comparison of CME- 
associated, coronal EUV dimmings observed simultaneously on the solar disk 
and above the solar limb. Both data sets yield a consistent picture of the eruption 
and mass loss during the early phase of the CME, and thus, support the view 
that above- the- limb dimmings and global-scale, on-disk dimmings essentially 
reflect the same process, namely, the density rarefication in the corona due to 
the expelled CME. 

We compared the dimming process with the kinematics of the CME and 
found that during the period of strongest CME acceleration (approx. 05:25- 
05:50 UT) the EUV intensity in both data sets dropped abruptly. Also, the 
mass loss rates, inferred from both on-disk and above-the-limb EUVI intensity 
observations, showed a distinct extremum in this period, suggesting that the bulk 
of the coronal mass was ejected, as the CME left the Sun. Moreover, the extrema 
of the CME acceleration profile and mass loss rates were only a few minutes 
apart in time, refiecting the close temporal relationship between the kinematics 
of the CME and the dimming process. Two hours after the eruption, both the 
EUVI 195 A mass loss rates and the Dopplergram derived from Fe xil 195.12 A 
observations with the Hinode/ElS instrument indicated that the loss of coronal 
plasma with a temperature of 1.5 MK had virtually stopped. 

The main characteristics of the dimmings in the event under study were 
comparable with other findings. For example, an intensity decrease within 20- 
30 min and a gradual recovery of the EUV brightness lasting several hours are 



characteristic features of the coronal dimming process {e.g., Aschwanden et al. 
19991 [Mcintosh et al.\ |2QQ7| [Reinard and Biesecker| [2^081 [Aschwanden et al. 
2009b'). Further, the observed amount of intensity drop, including the difference 
in brightness decrease between small-scale and global-scale dimmings, is typical 



for the two types of dimmings ( 


Aschwanden et al. 1999 Zarro et al. , 


1999 


Thompson et a/., 


2000 Mcintosh et al.^ 


2007 Harra et al.^ 


2007). Also, the 



initial growing of dimming areas as well as the apparent motion of these regions 



away from the source region of the eruption was reported earlier (Mcintosh et 
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a/., 2007). In the 13 February 2009 CME event, however, we observed for the 



first time an additional motion pattern of the small-scale twin dimming regions, 
namely, an apparent clockwise rotation around the source region of the eruption 
or the 'center' of the flare site, respectively. Taking the orientation of the first and 
last appearing postflare loops as a reference, this rotation can also be described 
as a 'sheared-to-potentiaP migration, superposed on the outward motion. 

How can we explain such an apparent rotation? Small-scale twin dimmings are 
often interpreted as the confined footpoint areas of a CME-associated, erupting 



flux rope (Sterling and Hudson 1997 Webb et al. 2000). The footpoints of a flux 



rope, however, are stationary. Therefore, moving bipolar dimmings represent not 
only the footpoints of the flux rope, but also other regions, where field lines are 
stretched. From MHD simulations the destabilization and subsequent eruption 
of flux ropes/CMEs is often achieved by twisting magnetic field structures (see 
e.g. |Chen et al\ |2002| [Fan and Gibson[ |2004[ [Torok and Kliem[ |2003[ ). Associ- 



ated dimming regions can as well be located at large distances from the source 
region of the eruption. For example, Delannee (2000) reported dimmings at the 
footpoints of transequatorial loops that connected a flaring active region to an 
active region located in the opposite hemisphere. They proposed that the appear- 
ance of such dimmings is strongly related to the magnetic field topology: When 
low- lying, active-region loops rise and erupt, thus forming a characteristic twin 
dimming pattern, the overlying large-scale loops rise and erupt as well.Mandrini] 
\et al. (^007) described a scenario in which the arcade above the erupting flux 
rope expands significantly before reconnecting. In this case, twin dimmings are 
expected to appear not only at the footpoints of the erupting flux rope, but also 
at the footpoints of the sheared magnetic arcade, as the expansion occurs. 

Bearing in mind this idea we propose a magnetic field topology in the event 
under study that may explain the observed rotation and outward migration of 
the small-scale twin dimming regions DRl and DR2 as well as the appearance of 



less and less sheared postflare loops, as the eruption progresses (c/. Figure 10). 
Before the eruption, a magnetic separatrix layer divides the flux rope from 
the overlying arcade that spans the polarity inversion line (PIL). The inner 
loops of this arcade, i.e., loops that are rooted closer to the PIL, are strongly 
sheared, low-lying loops, whereas the loop tops of the less and less sheared outer 
loops are at successively higher altitudes. Once the magnetic configuration has 
become unstable the flux rope erupts. As a consequence, the field lines of the 
flux rope are stretched and the first dimmings appear near its footpoints, ie., 
at strongly sheared positions. When the erupting flux rope encounters the field 
lines of the arcade, it stretches them while it continues to rise. Shortly before 
the stretched field lines reconnect, dimmings appear near their footpoints. Since 
the loop tops of the overlying arcade are located at different heights, the low- 
lying, inner loops, rooted at strongly sheared positions, are stretched before the 
large-scale, outer loops, whose footpoints are less and less sheared. Therefore, the 
respective dimmings appear at different times and positions, and the observer 
gets the impression of rotating twin dimming regions or dimmings that evolve 
from a 'sheared' to a more 'potential' configuration while they move away from 
the source region of the eruption, respectively. Once the stretched field lines of 
the overlying arcade reconnect, their upper parts supply the erupting flux rope 
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Figure 10. Apparent migration of the CoGs of the locahzed dimming regions. Top row: side 
view (a) and top view (b) before the eruption. The flux rope (hehcal structure) is rooted 
at both sides of the PIL (straight hue). Close to the PIL the overlying arcade consists of 
low-lying, strongly sheared loops, whereas at greater distances from the PIL the amount of 
shearing decreases and the altitude of the loop tops increases. Bottom row: side view (c) and 
top view (d) during the eruption. For clarity, only one fleld line of each fleld line triple from 
(a) and (b) is shown. The fleld lines of the overlying arcade are pushed upward by the erupting 
flux rope. Shortly before the thus stretched fleld lines reconnect, dim patches appear near their 
footpoints. The dark, fllled ellipses mark the positions of the CoGs of the coronal dimmings at 
different times (black/gray: CoGs near currently/recently stretched lines). The arrows point in 
the direction of the apparent motion of the CoGs. The white ellipses highlight the bright flare 
ribbons (energy deposition sites) at the footpoints of the currently reconnected fleld lines. 

with new poloidal flux, while their lower parts form the postflare loops. Since 
the first appearing postflare loops are created from the low- lying, inner, and 
strongly sheared loops, they are strongly sheared themselves, whereas the later 
appearing postflare loops are created from the less and less sheared outer loops, 
and thus, have a more and more potential configuration. 

In analyzing the behavior of the relatively small twin dimmings DRl and DR2 
on the solar surface we reach the conclusion that the observed motion pattern 
reflects not only the eruption of the flux rope, but also the ensuing stretching 
of the overlying arcade. The above- the- limb dimming region DR-L, however, is 
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larger than DRl and DR2 by a factor of 100 {cf. top panel of Figure [s] and 
Figure [sj^b)), and thus, rather represents the vast cross section of the coronal 
volume, evacuated by the CME. 

Finally, the existence of both STA and STB/EUVI observations of the 13 
February 2009 CME event allowed us for the first time to compare directly 
two methods of deriving the CME velocity in the initial stage of the eruption. 
The velocity values obtained from above-the-limb observations of the leading 
edge of the CME (CME tracking) and on-disk observations of the expanding 
global-scale dimmings were comparable (212 ± lOkms"^ and 228 ± 36kms~^, 
respectively). This demonstrates that global-scale, on-disk dimmings, originating 
from the source region of the eruption, propagate with a speed similar to that 
of the leaving CME front, and thus, supports the general view that global-scale 
EUV dimmings are on-disk signatures of CME launches. Accordingly, they can 
be used as a means for identifying CME source regions and estimating initial 
CME velocities, which is of particular interest in connection with Earth-directed 
halo CMEs and space weather forecast. 
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